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Abstract 

Sapqcd, the Abelian projection QCD(APQCD) action, is evalu- 
ated using the demon method. For SU(2), Sapqcd at strong cou- 
pling is essentially the compact QED(CQED) action with [icQED = 
\fisu(2)'^ extended and higher representation plaquettes are absent. 
Since CQED deconfines when [icQED >~ 1, this relation must break 
down as [isu{2) ~^ 2. Indeed we find Sapqcd mutates: near [isu{2) ~ 2 
it gains additional operators, including an exogenous negative mag- 
netic monopole mass shift. Sapqcd for SU{'i) has similar behavior. 
The Appendix gives a brief explanation of the demon method. 



A clear demonstration that monopole condensation is the origin of QCD 
confinement would be a notable achievement. To this end, 't Hooft [1, 2, 3] 
proposed that QCD monopoles are magnetic with respect to the 
Cartan subgroup of color SU{N). Full SU{N) gauge symmetry obscures 
these charges and it is necessary to gauge fix at least the SU{N) /[U{1)]^~^ 
symmetry to expose them. In this scenario monopoles are fixed-gauge man- 
ifestations of gauge field features responsible for QCD confinement. Only 
in special gauges does one have a picture of QCD confinement caused by 
monopole condensation. In other gauges the features causing confinement 
are still present but they do not look like magnetic monopoles [4]. 

Numerical studies have found that maximal Abelian(MA) gauge [5] is 
compelling for 't Hooft's hypothesis. Upon decomposing gauge field A into 
purely diagonal(n) and purely off-diagonal(c/i) parts 

A = A'' + (1) 

the MA gauge condition D'^Af = d^Af - igiA^.Af] = leaves a residual 
gauge invariance under 

N 

Residual = diag(exp-''-i , • • • , exp-''-~), ^ co, = 0. (2) 

8 = 1 

Under Oj-esiduai the N diagonal matrix elements (A")^^ transform as neutral 
photon fields whereas the N{N — 1) off diagonal matrix elements {A'^^)ij 
transform as charged matter fields: iA^)ii — )• {A^)ii — ^d^itOi and, for i ^ 
J-> {^t)^l {Af),jeyi^-'^^^-^i\ Since {A''^)^J carries two different U{1) 
charges, the A'^^ fields induce "interspecies" interactions between the N pho- 
tons. On the lattice the monopole currents are identified according to dis- 
cretized versions [6] of A;^ = j^ef^^xsd^fxs and /^j, = - d^A^^^. 

This procedure of using only the diagonal A^ components of the SU (N) 
gauge fields for measuring and f^^ is called Ahelian projection. Since 
trAj] = in SU{N), an irreducible representation of is 

^ i^h- (3) 
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Figure 1: To integrate out we: (i)generate an ensemble of impor- 

tance sampling SU{N) gauge configurations; (ii)project this ensemble to a 
[[/(f)]^~^ ensemble; and (iii)compute the Sapqcd couplings using the mi- 
crocanonical demon. 

The 9^ transform as — ^ — ^^t^'-^i ^i^d automatically obey constraint 
Ylf=i = 0. We shall refer to the quantum dynamics of the N angles 9^ 
as Abelian projected QCD or APQCD. Equivalently, APQCD is the field 
theory obtained by integrating out from QCD in MA gauge [7]. Its 
action Sapqcd is formally defined as 

- Sapqcd[9\ • • • , = log{ / [dA''] eM-SQCD) ^fp miK\] (4) 

where ^fp is the Faddeev-Popov determinant [8]. Monopoles arise in APQCD 
due to topological quantum fiuctuations in the compact fields 9\ 

While there is no guarantee that Sapqcd has a simple form or is other- 
wise well-behaved, it is of central import due to Abelian dominance [9], the 
fact that 9^ Wilson loops in APQCD have predominantly the same string 
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tension as SU{N) Wilson loops in QCD. Abelian dominance has the fol- 
lowing formal implication. If trW is an SU{N) Wilson loop and if {•)qcd 
and {•)apqcd refer respectively to Sqcd and Sapqcd expectation values, the 
APQCD operator W which obeys 

(W)apqcd = tr(W)QCD (5) 

is 

W = eM+SAPQCD) I exp{-SQCD) ^fp S[D;Af] trW. (6) 

Abelian dominance means that W, which in other gauges would be a com- 
plicated superposition of assorted [[/(l)]^~^-invariant operators of various 
sizes, shapes, and topologies, is (for string tension) well-approximated by a 
6^ loop of the same size and shape as trW in MA gauge. In other gauges the 
(W) APQCD string tension would be due to a combination of Sapqcd effects 
and details (such as operator coefficients) of W. In MA gauge, Sapqcd alone 
determines string tension: given Sapqcd one can reconstruct the QCD string 
tension using APQCD Wilson loops without reference to the RHS of (6). In 
this sense, in MA gauge Sapqcd knows the QCD string tension. 

Our numerical procedure for evaluating Sapqcd is schematically sum- 
marized in Figure I. First, focusing temporarily on SU{2) we make an en- 
semble of importance sampling APQCD gauge configurations by applying 
the Abelian projection to a set of Monte Carlo SU{2) gauge configurations 
at some coupling fisu(2)- Then, seeking the action Sapqcd which would re- 
produce this APQCD ensemble [7] in a Monte Carlo simulation, we state an 
ansatz for Sapqcd and apply the microcanonical demon technique [10] to 
compute the parameters of this ansatz. This whole "inverse Monte Carlo" 
procedure is repeated starting from different QCD ensembles to determine 
how Sapqcd varies with I3su(2)- 

The general [/(I)-invariant action consistent with APQCD symmetries 
involves an infinity of operators. However, previous studies [7, II] and, in- 
dependently, the demon technique indicate that neither extended nor highly 



3 



Table 1: S^^'^^^ couplings in APQCD measured by the demon on a 14"^ X 4 
lattice, where QCD is confining over the whole fisu(2) range given. 



parameter 




L = 1 


L -- 


= 2 


L -- 


= 3 




1.0 


.50(.02) 


mi 


.01) 


.02( 


.01) 


HL) 




-.02(.01) 


.00( 


.00) 


.01( 


.01) 


HL) 




.OO(.OO) 


.00( 


.00) 


.00( 


.01) 


fh{L) 


2.0 


1.03(.04) 


.00( 


.00) 


.02( 


.01) 


HL) 




.02(.00) 


.00( 


.00) 


.00( 


.01) 


HL) 




-.Ol(.Ol) 


-.OK 


:.oi) 


.00( 


.00) 


HL) 


2.2 


1.2(.03) 


.00( 


.00) 


.02( 


.01) 


HL) 




.04(.01) 


.00( 


.00) 


.00( 


.00) 


HL) 




-.Ol(.Ol) 


-.03( 


:.o2) 


-.02( 


:.oi) 



charged Wilson loops contribute substantially to Sapqcd- In particular, Ta- 
ble 1 shows the results of applying the demon to the ansatz 

- ST"-'^ ^jzjZHL)T. cos q<d,.{L) (7) 

L=l g = l x,fl<v 

where cos qQf^i,{L) is an L X L plaquette in U{1) representation given in 
terms of link angles fiq{L) is its coupling the demon computes. Appendix 
A explains the technical idea underlying the demon technique. Figuratively, 
imagine a battalion of demons each carrying M coupled thermometers whose 
temperatures correspond to the M undetermined couplings in the ansatz 
action. {M = 9 for S'^''^^^.) The demons thermahze with an APQCD 
configuration — the heat bath — by hopping from link to link. Each hopped- 
upon link is randomly updated, and the hopping demon absorbs or emits 
the corresponding energy(=action) decrement or increment created in the 
configuration by the update. The thermometers are coupled by requiring all 
their energies to simultaneously remain within a given range [— if 
any proposed update violates this range it is rejected. Upon thermalization 
the couplings are read off by fitting the energy distributions of the demons' 
thermometers to M Boltzmann distributions. Statistical errors are computed 
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by jackknifing the demons. (The errors from jackknifing SU{2) configura- 
tions are comparable.) In principle, if S^^'^^^ contains all the operators of 
SapqcD} the demon method tells us which couplings vanish and measures 
all the nonzero couplings exactly (modulo statistics). In practice, S^^'^^^ is a 
truncated action which is unlikely to contain all Sapqcd operators. Extensive 
numerical experiments with control ensembles (see Appendix and Ref. [15]) 
reveal that if important operators are missing the method yields readjusted 
"effective" coupling values. These effective values are not equal to the true 
values. 

As illustrated in Table 1, L = q = 1 plaquettes dominate Sapqcd- 
Furthermore, there is no significant signal for any of the L > 1 plaquette 
couplings at any fisu(2)- This result, unanticipated, has a substantial impli- 
cation. Since BKT transformations of Z = 1 actions do not lead to extended 
monopoles, this implies that, at least within our fisu(2) range, the fun- 
damental dynamical degrees of freedom [11] for confinement in APQCD are 
pointlike P rather than extended monopoles, as considered in [12]. In 
future work, it will be important to examine if this result survives the zero 
lattice spacing limit, e.g. if APQCD monopoles truly are pointlike. 

On the other hand, the L = 1, q = 2 plaquette has a noticeable signal 
at fisu(2) = 2.2 in Table 1. Therefore, we focus now on the L = 1 ansatz^ 

- ST't- = E E /3,cos(?0^. - Kj2k,{x)k,{x). (8) 

q = l :r,/i<i^ x,i^ 

The K, operator shifts the q = 1^ 1^ monopole mass [12] implicit in /3i cos O^j,, 
allowing the APQCD monopole mass to be independent of Of course, 
l3q and k, vary with fisu(2)- Figure 2 shows S^^^^"^ coefficients [^2 and /€, 
computed by the demon, as a function of fisu(2)- l/^sl, not depicted, is always 
smaller than |/32|, typically by a factor of 3 — 5. Fach fisu(2) configuration 
is generated fresh from a cold start so our data points do not contain any 
spurious correlations. Our Ng X Nj = 20"^ X 16 lattices are all well inside the 

"'^ Unless otherwise specified, _L = 1 is assumed in the remainder of this Note. 
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1 2 

/^SU(2) 

Figure 2: Figure 2 depicts 5*^"^^*^ coefficients [^2 and function of 

Psu(2)- l/^sl, not depicted, is always smaller than |/32|, typically by a factor 
of 3 — 5. Our 20"^ x 16 lattices are all well inside the zero temperature phase 
for the [isu(2) range depicted. The bold /3i = \[isu(2) line is a guide-to-eye. 
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zero temperature phase; for the range of fisu(2) shown the APQCD Polyakov 
loop vanishes. As depicted, at strong couphng(/35(7(2) < 2) 

/3l ~ ^/35C/(2), /32,3 ~ 0, /€ ~ 0, (9) 

that is, Sapqcd reduces to the compact QED(CQED) action at strong cou- 
phng. At weaker couphng(/35[/(2) > 2) [^2 and n grow in magnitude but /3i 
always remains the largest coupling. 

Note that since monopoles are condensed when ficQED < 1 in CQED [16], 
Figure 2 or Eq. (9) vicariously proves that SU (2) monopoles are condensed 
when I3su(2) < 2. 

When fisu(2) > 2, the situation is not so clear. In fact, Eigure 2 suggests 
a paradox in the fisu(2) > 2 region: how can APQCD maintain confinement 
in the continuum limit if CQED deconfines when ficQED > 1? Clearly, 
either the meaning or validity of relation (9) must break down when fisu(2) is 
sufficiently large. Either (I)Abelian dominance does not survive the fisu(2) ~ 
2 crossover making Sapqcd less pertinent at weaker coupling — see discussion 
pertaining to Eq. (6); or {1\-)Sapqcd gains additional operators near fisu(2) ~ 
2; or a combination of (I) and (II). We do not have anything to say about (I) 
in this Note except that Abelian dominance apparently has been observed 
at all ^su(2) < 2.6 [9]. 

(II) requires that when fisu(2) > 2 Sapqcd is no longer well described 
by ScQED- Indeed, we can demonstrate this by simulating 

- ScQED = X] I^CQEDCOsQ^i, (10) 

^'^<'' \I3CQED=MPSU(2)) 

(also on a 20"^ X 16 lattice) to see if it reproduces corresponding APQCD 
expectation values. As depicted in Eigure 3, Scqed reproduces APQCD pla- 
quette averages and monopole densities only in the SU (2) strong coupling 
region. At weaker coupling the CQED simulations start to disagree dra- 
matically with APQCD. This implies that at weaker coupling either other 
terms of S'^^^^'^ have become important or S'^^^^'^ itself is inadequate. In 
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Figure 3: 3A compares APQCD plaquettes Papqcd at fisu(2) to CQED pla- 
quettes Pcqed at /3cqed = /3i{/3su(2)) for a range of /3su(2) values. When 
Psu(2) < 2 the data points he on the bold Pcqed = Papqcd line showing 
that ScQED is a good model of Sapqcd- The set of points lying off of the 
Pcqed = Papqcd hue corresponds to /3su(2) > 2, when Scqed is not a good 
model of Sapqcd- 3B is an analogous plot using monopole densities. 
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any case, this means Sapqcd is not form-invariant between the strong and 
weak couphng regimes: at strong couphng Sapqcd is well approximated by 
ScQED] at crossover region fisu(2) ~ 2 Sapqcd mutates and develops sub- 
stantial deviations from Scqed- Inspection of Figure 2 reveals that a possible 
scenario^ might be that /€, the exogenous magnetic monopole mass shift, be- 
comes more and more negative at larger fisu(2)- As negative monopole mass 
favors monopole condensation (compensating for a large the occurrence 
of a sufficiently negative n in Sapqcd at fisu(2) > 2 could maintain APQCD 
confinement. 

Note that Figure 2, as characterized by Fq. (9), "explains" Abelian 
dominance — at least in the strong coupling regime. The SU (2) plaquette in 
the strong coupling expansion behaves like Pqcd ~ \[^su(2) and the CQFD 
plaquette hke PcQSD ~ \I^cqed- Therefore, identifying PcqsdI/^cqsd = A) 
with Papqcd and applying Fq. (9) yields 

Papqcd ~ -^f^su(2) ~ Pqcd- (H) 

Carrying this argument over to larger Wilson loops leads to a strong coupling 
version"^ of Abelian dominance: at sufficiently strong coupling APQCD and 
QCD Wilson loop averages and, hence, string tensions are equal. Figure 4 
confirms (II) and shows how this relation breaks down at weaker coupling. 
Interestingly, Fq. (II) contradicts the naive expectation, based on Pqcd 
containing a trace over a 2 X 2 matrix and Papqcd involving no trace, that 

Papqcd = \Pqcd- 

We have obtained similar results for the SU (3) Abelian projection which 
will be described elsewhere. For S'[/(3), Sapqcd is more complicated due to 
interspecies dynamics [II, 17, 18]. Nonetheless, we have observed completely 
analogous behavior in SU{'i). At stronger couplings, Sapqcd is dominated 
hy L = q = 1 operators; at weaker couplings, there is a crossover to a more 

^Our demon studies, exemplified in Table 1, seem to rule out the alternative possibility 
that L > 1 plaquettes become important at larger couplings. 

^Not to be confused with weak coupling Abelian dominance which requires only string 
tension equality. 
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Figure 4: Figure 4 depicts the APQCD and SU{2) plaquett function 
of fisu{2)- In the strong couphng vegion{fi su {2) < 2), both the APQCD and 
SU{2) plaquettes grow hke j[Ssu{2)^ the guide-to-eye hue's slope. At weaker 
couphng(/35(7(2) > 2) the APQCD plaquettes deviate substantially from the 
SU{2) plaquettes. Correspondingly, the monopole density decelerates no- 
ticeably near fisu(2) ~ 2. 



complicated, but still L = 1 action. Again, L > 1 plaquettes are never 
resolved. Very preliminary SU{3) results are reported in Ref. [3]. 

In conclusion, while it is not form-invariant between the strong and 
weak coupling regimes, Sapqcd is dominated by Z = 1 operators for all 
values of fisu(2) studied. This implies that P (rather than L'^) monopoles 
are the dominant fundamental dynamical degrees of freedom for APQCD 
confinement, and that phenomenological features of the APQCD confinement 
mechanism, such as whether APQCD is a Type I or Type II superconductor, 
might vary with lattice spacing. In particular, the results of this paper predict 
APQCD is a Type II superconductor like CQED [19] when fisu(2) is small, 
as distinguished from the fisu(2) > 2 case studied in [20]. 

4. Acknowledgements 

It is a pleasure to thank Dick Haymaker and Misha Polikarpov for stimu- 
lating discussions and MP for his SU{2) FORTRAN codes. I am indebted to 
the Institute for Theoretical and Experimental Physics(ITEP) for their hos- 
pitality. Computing was done at the NERSC Supercomputer Center. The 
author is supported by DOE grant DE-EG05-91ER40617. 

References 

[1] G. 't Hooft, Nucl. Phys. B190 (1981) 455. 

[2] A. Kronfeld, G. Schierholz, U. Wiese, Nucl. Phys. B293 (1987) 461. 

[3] More references are given in K. Yee, Helv. Phys. Acta 67 (1994) 345. 

[4] An illustration is suggested in M.N. Chernodub, M.I. Polikarpov, A.I. 
Veselov, preprint ITEP-63-94 (Aug., 1994). M.I. Polikarpov, private 
communication. 

[5] A. Kronfeld, M. Laursen, G. Schierholz, U. Wiese, Phys. Lett. B198 
(1987) 516. 



11 



[6] T. Banks, R. Myerson and J. Kogut, Nucl. Phys. B129 (1977) 493; D. 
Toussaint and T. DeGrand, Phys. Rev. D22 (1980) 2478. 

[7] K. Yee, Phys. Rev. D49 (1994) 2574. 

[8] The Faddeev-Popov determinant is discussed in Ref. [4] and M. Pohkar- 
pov and K. Yee, Phys. Lett. B333 (1994) 452. 

[9] T. Suzuki and I. Yotsuyanagi, Phys. Rev. D42 (1990) 4257; S. Hioki, S. 
Kitahara, S. Kiura, Y. Matsubara, 0. Miyamura, S. Ohno, T. Suzuki, 
Phys. Lett. B272 (1991) 326. 

[10] M. Creutz, A. Gocksch, M. Ogilvie, and M. Okawa, Phys. Rev. Lett. 53 
(1984) 875; M. Creutz, Phys. Rev. Lett. 50 (1983) 1411. 

[11] K. Yee, Phys. Rev. D50 (1994) 2309. 

[12] J. Smit and A. van der Sijs, NucL Phys. B355 (1991) 603. Refs. [11, 13, 
14] give some monopole and extended monopole densities. 

[13] V. Bornyakov, E. Ilgenfritz, M. Laursen, V. Mitrjushkin, M. MuUer- 
Preussker, A. van der Sijs, A. Zadorozhny, Phys. Lett. B261 (1991) 116. 

[14] G.L Pouhs, H.D. Trottier, R.M. Woloshyn, TRIUMF preprint TRLPP- 
93-88 (revised, 1994). 

[15] K. Yee, manuscript in preparation. 

[16] M.L Pohkarpov, L. PoUey, U.J. Wiese, Phys. Lett. B253 (1991) 212. 
[17] M. L Pohkarpov and K. Yee, Phys. Lett. B316 (1993) 333. 
[18] K. Yee, Mod. Phys. Lett. A9 (1994) 1991. 

[19] V. Singh, D. Browne, R. Haymaker, Phys. Rev. D47 (1993) 1715. 
[20] V. Singh, D. Browne, R. Haymaker, Phys. Lett. B306 (1993) 115. 



12 



Appendix A The Demon Method 



This Appendix explains the demon method [10] lor ansatzes which can be 
link- wise partitioned so that S*^""^*" = E^,^ ^^^^^ ^^7^^^ = ^'^^ 

has the same form and coupling values jSi on all links, i labels the M different 
operators Hi. While all the operators of .S*^^^^*^ share a common link Ux,^^ 
they may (and do) depend on other links. For S^^'^^^ of Eq. (8), Af = 4; 
Hi, H2, and H3 correspond to oc Y^i,=-a cos qQy,i,{x) for (j' = 1, 2, 3; and H^^ is 
oc Xlj/ k^[x)k^[x) summed over all directions affected by link 6j^. Let 

n = n{Hi,H2,---) (A.i) 

denote the number density of states in the "energy" interval Hi and Hi + Si. 
In statistical mechanics language, the entropy is proportional to log and, 
for a system with fixed total energies Ef , the inverse temperatures are 

A - ^ • (A.2) 



OH 



T 



Now imagine a "demon" carrying M thermometers corresponding to the 
M operators Hi. The job of the demon is to measure the inverse temperatures 
l3i of a heat bath — in our case an APQCD gauge configuration. To do this, 
the demon hops link-to-link and exchanges energy with the bath until it ther- 
malizes. More precisely, at each link the demon thermometers are changed by 
energy increment computed as follows: the bath link is randomly 

updated and Ai7f^™°" = Hf'^ - is computed by evaluating Hf'^ and 

on the bath before and after the update. Since Ai^f^*^ = _Ai7f^™°", 
the demon plus bath energy is constant under this procedure.'* The ther- 
mometers are coupled by requiring every thermometer energy to be inside 



■*In practice we do not retain the update of the APQCD configurations, so that the 
demon plus bath energy is not really constant in our procedure (as it is in Ref. [10]). 
Nonretention shortens the computer algorithm and avoids any possibility of damaging the 
APQCD configuration, a real danger since we have a whole battalion of energy-absorbing 
demons. 
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some range [— -E^, -E^]; if a potential update pushes any thermometer outside 
this range, it is rejected. 

Upon denoting heat bath quantities with primes, the microcanonical 
partition function of the total demon-bath system is 



zueLe^,---) = j[de][de'] f{ 8{h, + hi - Et) 

A=l 

. M 

= J[d9][dH']n{H[,H^,---) HSiH^ + Hl-Ef) 



A=l 

Performing the [dH'] integration and Taylor expanding entropy log yields 

,. M 

ZUEL E^, ■■■) = ^{E^, E^, ■■■) [dO] exp ^ -[i,H,. (A.3) 

A=l 

Eq. (A.3) expresses the well known result that a thermalized subsystem of 
a microcanonical ensemble has a Boltzmann distribution with inverse tem- 
peratures given by (A. 2). Hence, if we set a battalion of demons free in 
the importance sampling APQCD configurations, upon thermalization the 
demons will return with their M thermometers each distributed in a Boltz- 
mann distribution. Therefore, the .S*^^^^*^ coupling constants fii corresponding 
to APQCD are readily extracted by fitting each of these M distributions to 
exp{-/3ji7j. 

We have tested the demon method extensively as follows. First, we gen- 
erate an ensemble of U{1) configurations according to a known U{1) action 
5*0. For example, 5*0 m&iy be ^^^^^^^ 

some point P in fiq{E) space. Then we 
apply the demon method with a trial ansatz Sq^^^^"^ which, for purposes of 
discussion, may or may not contain all operators of true action 5*0. If Sq^^^^"^ 
contains all operators of 5*0, then we find that the demon always successfully 
recovers P and 5*0, that is: (i)coefficients of operators in Sq^^^^"^ which do 
not exist in 5*0 vanish modulo statistical errors comparable in size to those 
in Table 1; and (ii)coefficients of operators in 5*0 equal P modulo statistical 
errors. The ability of the demon method to unambiguously and automatically 
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reveal when an operator does not exist in 5*0 is an advantage ol the demon 
method over other methods. 

II does not contain all the operators of 5*0 , the situation is less 

straightforward. The demon apparently tries to obtain an effective action by 
using the available operators in to fit the ensemble as optimally as 

possible. However, is not easy to nail down what exactly is being optimized. 
Therefore, it is important to simulate the ansatz action with demon couplings 
and verify, as we have, that ansatz expectation values reproduce APQCD 
expectation values. 
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